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I. INTRODUCTION 



The Be is the only heavy meson consists of two heavy quarks with different flavors, hence 
decay properties of this meson are of special interest. The difference of heavy quark flavors 
forbids annihilation of this meson into gluons, as a result of which the excited Be meson states 
lying below the the threshold of decay into the pair of heavy B and D mesons undergo pionic 
or radiative transitions to the pseudoscalar ground state which is more stable compared to 
the corresponding quarkonia and decays mostly weakly. It is expected that the experimental 
studies on the Be meson and its decay properties constitute one of the important parts 
of physics program at LHC The study of heavy mesons will provide not only a window 
in extracting the most accurate values of Cabbibo- Kobayashi- Maskawa(CKM) matrix 
elements as the sources of CP-violation in Standard Model (SM) but also in understanding 
the perturbative and non-perturbative aspects of QCD. 

In the present work, we work out the rare radiative Be — ?■ Dgi (2460)7 transition in 
the framework of QCD sum rules [l, 2|. Here, the Dsi{2A60) is the axial vector charmed- 
strange meson with quantum numbers = 1^ and interpolating current r/j, = 37^750. 
This transition proceeds via both weak annihilation (WA) and electromagnetic penguin 
(EP) of flavor changing neutral current (FCNC) transition, based on 6 — > 57 at quark 
level. Using the quark condensate, mixed and two-gluon condensate diagrams as well as 
propagation of the soft quark in the electromagnetic field as non-perturbative effects, we 
calculate the transition form factors responsible for this decay in both weak annihilation 
and electromagnetic penguin modes. We then use these form factors to estimate the decay 
rates and branching ratios in these modes as well as the total decay rate and branching 
fraction. As it is expected, the dominant contribution comes from the weak annihilation 
part. Note that the Be — ?■ -D*7 transition is studied in the same framework in [3]. Some 
other radiative channels of the Be meson like Be IV'j and Be — > i?*7 have also previously 
studied using QCD sum rules technique {4, 5|. For some other decay channels of the Be 
meson see for instance 

The outline of the paper is as follows. In next section, considering the radiation of the 



photon from both Be and Dsi mesons, we construct the transition amplitude for the weak 

( B ) 

annihilation channel in terms of four relevant form factors. Two of the form factors Fy " 

( R ) 

and F\ " responsible for the emission of the photon from the initial state are calculated 
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FIG. 1. The weak annihilation mechanism for Br 



in jj]. The remaining two form factors, Fy ^ and " representing the emission of the 
photon from Dgi meson are calculated in section III . In section IV, considering the two 
gluon condensate contributions, we calculate the transition form factors responsible for the 
electromagnetic penguin mode. Finally, section V is devoted to the numerical analysis of 
the form factors, calculation of the decay rates and branching ratios for the modes under 
consideration together with total decay rate and branching ratio. This section also contains 
our concluding remarks. 



II. WEAK ANNIHILATION AMPLITUDE 

In this section, we construct the weak annihilation amplitude for the radiative Dsil 
transition. Considering the quark contents of the initial and final mesonic states, the possible 
diagrams are shown in figure 1. Taking into account these diagrams, the transition amplitude 
for the radiative decay under consideration is written as 

M^^(i?, ^ Z},i7) = ^K.V;:(D,i(p)7(g)|(^r,c)(cF'^6)|5,(p + g)), (1) 

where Gp is the Fermi weak coupling constant, Vij are elements of the CKM matrix, F,^ = 
7^(1 — 75); and p, q and p+q are the momenta of the Dgi, photon and Be meson, respectively. 
To proceed further, we use the factorization hypothesis and write the matrix elements in 
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Eq.p as 



where we have divided the matrix elements into two separate parts: emission of the photon 
from Be meson (diagrams i and ii in figure 1), represented by the covariant decomposition 
Tj},^"^ and the emission of the photon from Dsi meson denoted by the decomposition 



(Dsi) 



(see diagrams iii and iv in figure 1). In Eq.([2]), the /^^ and fn^^ are the decay constants of the 
Be and Dsi mesons, respectively. The and e^^"'^^'^ stand respectively for the polarization 
vectors of the photon and Dgi meson. The covariant decompositions T^f"'' and Tjiu"^^ are 
defined in the following way: 



T!,^^\p,q) = ^ I rf'xe^^^X0|r|j7(x)c(0)r,6(0)U,(p + g)) 



(3) 



and 



(4) 



where j^™ is the electromagnetic current and T is the time ordering operator. Applying the 
Ward identity for the electromagnetic current and using = for the real photon as well as 
e.q = and e^^^^Kp = similar to what is done in [3| and [l^, we get the following results 
correspond to the emission of the photon from the initial and final mesonic states in terms 
of form factors: 



(£.£(^»i))(p.g) - {e.p) [e^^^'lq 



iF 



(Be) 



(5) 



iF 



{Dsi) 



J- f m (c ciDsl) \ _|_ e- AD al) V cptj>^„'y rpiDsl) I 

+ jDsiT^Dsi K-^ ) + £u^l\a^ B p q ry > 



(6) 



where F^'^^ and F^^? are the transition form factors. Using these expressions, finally, we 



V(A) - V(A) 

get the WA transition amplitude as 



Be 



(7) 



As we previously mentioned, the form factors Fy " and F\ " corresponding to the emission 
of the photon from Be meson are calculated in ^. We calculate the remaining two form 



factors Fy^"^'^ and F^"^' in the next section. 



III. QCD SUM RULES FOR THE FORM FACTORS F^^''^ AND F'f'^^ 



To calculate the transition form factors in QCD sum rules formalism, we start considering 
the following correlation function: 

U^Ap,q)=^ j rfW«^(7(g)|r|c(a;)7^(l-75)s(x)s(0)7.75c(0)}|0), (8) 

where Q = p + q. The basic idea in this method is to calculate this correlation function once 
in hadronic language called phenomenological or physical side and then in terms of QCD 
degrees of freedom using the operator product expansion in deep Euclidean space, called 
the theoretical or QCD side. The two representations are then matched in order to get the 
QCD sum rules for the form factors. To suppress the contributions coming from the higher 
states and continuum, we apply Borel transformation as well as continuum subtraction. 
These procedures bring two auxiliary parameters, namely Borel mass parameter as well 
as continuum threshold that we shall find their working regions such that the physical 
observables be independent of these parameters in these regions. 

First we focus on calculation of the phenomenological side. For this aim, we insert a full 
set of hadronic Dgi state to Eq.([H]). After performing integral over d'^x, we get 

{l{q)\ci^{l - l5)s\Dsi{p)){Dsi{p)\s-f^-r5c\0) 



^t^u{p,q) 



(9) 



The matrix element, (-Dsi(p)|s7;y75c|0) is defined in terms of decay constant and polarization 
vector of Dsi meson, 



(10) 



while the transition matrix element is parametrized in terms of form factors. 



{l{<l)\cl^l{^-lb)s\Dsl{p)) 



m 



+ 



m 




{Dsi) 



Putting all above equations together and summing over the polarization vector of the D^i 
meson, we get the following result for the phenomenological part of the correlation function: 



+ 



The QCD side of the correlation function is calculated in deep Euclidean region in terms 
of QCD parameters. To proceed, first we write the correlation function in terms of two 
selected structures as 



m, 



m 



[12] 



(13) 



where each function Ilj = 1 or 2) has perturbative and non-perturbative parts, i.e., 

n,; = nf + nr""^""*. (i4) 

For calculation of the perturbative parts, we consider (a) and (b) in figure 2 where photon 
can be radiated from both charm and strange quarks. For the non-perturbative effects, we 
take into account the quark condensate and mixed diagrams [Figs. 2(c), 2(d) and 2(e)] as 
well as diagram 2(f) indicating the interaction of the photon with the soft quark. 
The perturbative part in each case can be written via the dispersion relation as 

/P'is 
ds— — + subtraction terms. (15) 
s - Q2 

where pi are the spectral densities. Our main task in this section is to calculate these 
spectral densities using the diagrams (a) and (b) in figure 2. Here we use a Feynman and 
Shwinger parameterization based method with several Borel transformations (see also 11|). 
The Feynman amplitude for the diagram (a) can be written as 



+ ml 



-7/.(l - 75) 



{Q + ky — rnl {p + ky — ml 
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FIG. 2. Diagrams for bare-loop [(a), (b)], power corrections from the operators with 3, 4 and 5 
dimensions [(c), (d), (e)] and propagation of the soft quark in electromagnetic field (f). 

Using first the Feynman parameterization, we perform integral over d'^k. Then, we use the 
Shwinger parameterization as 

1 1 '•^ 



daa"' e 



n—l„—aA 



(17) 



A" r(n) 7o 

to write the denominators of the obtained results in exponential forms. As a result, we get 



■pert 



T-fper 



1 /■! 



^^^y^-j I dxx J dy ms{'mc + msXy)+p'^xx{l — xy) + 2p.qxx^y^ I dae 



If 



-aA 



Yfpert gNcQ^ 



dyx ms{mc — msxy) — xx{l — xy) — 2p.qxx^y^ 



f 

Jo 



(18) 
(19) 



where x{y) — 1 — x{y), and A = m^x + m^x — p xx y — Q xxy. 
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Applying double Borel transformation B{Mf)B{M^) on nf *, that transforms Mf 
and — > M|, we obtain 

ffpert _ eNcQs OXOI C A ^ (mgx+mgx)(ai+a2) 



47r2 (Ji + (72 Jo X 

fn(7i 2— /''^l Cr2\ _2'^? I 

rricms + m x h p xa; I a; 1 ) + zp.qxx — > , 
Cri+Cr2 V (7i + (72 (7i + (72 / ((7i + (72)"' 



(20) 



- 47r2 a, + (72 X 



2 CTl 2- "^l , 0-2 A „ _ 2 CTl 



^ — . ^ , - ) - 2p.qxx'^ . , 

(7i + (72 V (Ti + (T2 (Ti + (72 / ((7i + (72j 

'(21) 

where (7i,2 = I/MI2 and we have used the Borel transformation of exponential function as 

Bp2{M^)e-''P' = 5{1 - aM^). (22) 

Now, we perform the second double Borle transformation on TL^^^* in order to transform ai 
and (72 to new variables s and t, using 

1-/1 \ - /I \ fr^'"'* 

ft(s,t,/) = B(-,aAB(-,a2)^, (23) 



St \S / \t / (7i(72 

where 

S(i (71)5(1,(72) e-^^'^^^-^) = (5(1 - ^(1 - j), (24) 



and we have also used 



^n^-aa ^ (_^)ng-aa_ (25) 

da 



The final expressions for the spectral densities are calculated via the following formula: 

Pi[s,P = / dt— —. (26) 

J t-p^ 

After lengthy calculations, we get the following spectral densities corresponding to the dia- 



gram (a): 
Pia(s,P^) 



dx—l mcms{s — p^)^x — 2ml{ml — p^)x^ 



47r2 (s-p^y,,^^ 

Am^x + p^ ^4p^ + — p^)^xx + ml ^(s — p^) + 4p^a;^ 
- 2m^x + (^s^ + p^s{x -l)+p^{x + 2)^ 



+mlx 



+x 



-m 



+ x^) [{s - p^) + 2p^x^ 



(4x2 - 2x)^ 



(27) 



P2a[S,P 



eNcQs 



An^ (s - 

+x 



— ^ / dx—2<mcms{s-p'^)'^x + 2m^{rril-p^)x'^ 



2mlx - (^s^ + p^s{x - 2) + ^"^(0; + 1)) 
+m^p2(-l + x^) (^(s - + 2p2x) -ml({s- p^) + ^^(4x2 - 2a;)) 

—m^x — 4mlx + p^{3p^ + s)xx + ml (^{s — p^) — Ap^x^^ 

where the integral boundaries xq and xi satisfy the the following inequality: 



(28) 



sxx — {mix + mix) > 0, 



(29) 



that comes from the Heaviside-Theta function arising in calculations. Similarly, we calculate 
the contribution of diagram 2 (b) . The final expressions for the spectral densities correspond 
to the two selected structures are obtained as 

eiV, 



Pi{s,p ) 



247r2 (p2"37)l|^^ (gsMI -2a- 2{a - ^f) + Sp^s^fl - Aa^ - 2a{?> + ^)(1 + 2/3) 
+/3(2 + 5/3) + a[2 + a(17 + 8/3)]} + p^(5 + 2a'' + a(5 - 4/3) - 7/3 + 2/3") 
-pS[4 + Aa^ -5p- 11/3^ - a^(19 + 8/3) + a(23 + 26/3 + 4^^)]) -6(mcms(p^ - s)^ 

+sa[2p^ + s'{-l + Aa){a - /3) +^^^(2 -7a + 3/3)]) In ( ||^~|^~^ )" 

+gc X(6s^{l - 2(a - p)^ - 2p)p + 3p2s2(l + a{2 + 5a) - 2(3 + a)(l + 2a)p - Ap^ 

+/3(2 + (17 + 8«)/3)) + p6(5 + 2/32 + /3(5 - 4a) - 7a + 2a') - p^s[4 + 4/3^ 
-So; - 11q;2 - I3'{19 + 8a) + /3(23 + 26a + Aa')]J - 6(mcms{p' - s)' 

+sl3[2p' + s'{-a + /3)(-l + 4/3) + p's{2 + 3a- 7/3)]) In (1^^±|^)] \, (30) 



and 



P2(S,P^) 



eNr. 



247r2 (p2 _ s)3 1^^^ [^(^5^5(4 + a(23 + a(-19 + 4a)) - 5/3 + 2(13 - 4a)a/3 
+(-11 + 4a)/32) + 3pV(-l + (-4 + a)a{-l + 4q;) - 2/3 + 2(7 - 4a)Q;^ 
+ (-5 + 4a)/32) - p6(5 + 2a^ + /3(5 + 2/3) - a{7 + 4/3)) 
+6s^q;[-1 + 2{a + - 2q;^ + /3^)]) - 6(|mcm,(p^ - s)^ 
-sa[2p^ + s'(-l + 4a)(a - /3) + ^^^(2 - 7a + 3;g)]) in + ~ ^ ~ 



+Q 



l + a-p + x). 

A (^6s^^(-l + 2(a - ^)2 + 2p) - p%5 + a(5 + 2a) -7p- Aap + 2p^) 
+pS{4 + a2(-ll + 4/3) + a(-5 + 26/3 - 8/3^) + /3[23 + /3(-19 + 4/3)]} 
+3p2s2[-l + a\-5 + 4/3) + (-4 + /3)/3(-l + 4/3) - 2a(l - 7/3 + 4/3')]) 
-6(^mcm,(p^ - s)^ + sp[-2p^ + s^(a - /3)(-l + 4/3) 

V.(-2-3..7«l)l.(i^^)]}. (31) 



where a = ^, ^ = ^ and A = yrTa2T^9232^r^2;9^^2^. 

For the non-perturbative part, we first calculate contributions of the quark condensate 
and mixed diagrams ( diagrams c,d and e in figure 2). After lengthy calculations, we obtain 
the result as follows: 



Yjnon—pert 

l(c,d,e) ~ 



7mj. 



Ami 



.2^2 ^4^2 ^2^4 ^4^4 



+ - 



{ss) 



2ml 



8ml 



+ 



2ml 



3mr 



+ 



2ml 



mr 



"^0 /- \ 



r2i?6 r^Ri ' r^R4 r'^R^ ' r^R-^ r^Ri\ 

3mc 1 



-6ml 24ml QitlI 



+ 



imr 



6ml 



Ami 



^2^6 ^6R4 ^4^4 ^2^4 ^6^2 ^4^2 J ^2^4 



(32) 



and 



T-rnon—pert 
^^2(c,d,e) 



rrir. 



,2 



{ss) + —{ss) 



2 



+ ^{ss) 



"^0 /- \ 



2 
2ml 

- 

2ml 



1 



1 



Am . 



3mz 



m„ 



i?4 r^Ri r^Ri 2r4i?2j 



+ 



8ml 



2ml 



+ 



3mr 



2ml 



+ 



mr 



im 



r^R'^ r^R'^ r'^R^ r^R^ r^R^ 
^ 2ml Amr 2ml 



mr 



^2^6 ^6^4 ^4^4 ^2^4 ^6^2 ^4^2 



+ 



Ami 



{ss), 



(33) 



where — — mi and B? — Q"^ — 



ni,. 
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The final contribution to the WA mode is that of diagram (f). This diagram corresponds 
to the propagation of the soft quark in the external electromagnetic field. Here we need to 
use the light-cone version of the QCD sum rules and photon distribution amplitudes (DAs). 
The relevant correlation function is of the form: 

n^.,{/)(p,g) =^ / f/'xe-^«^(7(g)|T{s(0)7;.75c(0)c(x)7.(l-75)s(x)}|0). (34) 



Contracting c-quark lines in Eq. fH9|) and using the propagator of the heavy quark in mo- 
mentum space, we obtain: 

^^^ (27r)4 m? - A;2 + m^^ui^ - 75)s|0). (35) 

To relate the matrix element exists in the above equation to the photon DAs, we use the 
identities, 

i 

Ifilalu = gfialu + gualf, " gfiula + iS f^vcXlXlb- (36) 



The relevant photon DAs of twist 2, 3, and 4 12, ll3|] are as follows 



(7(?)k7;.'S|0) = -%/37 I dw/j^^\u)x'^ Fe^{ux) 



2 











{l{q)\s(Tai3s\0) = Qs{ss) / du(j){u)Fap{ux) 

Jo 



16 Jo 

^QJ^ f duB{u)x''{xpF^p{ux)-x^Fpp), (37) 
^ Jo 

where F^^ is the field strength tensor of the electromagnetic field and is defined by 

F^u{x) = ]^e^yapFaf}{x). (38) 

The wave function 0(m) is defined in terms of the magnetic susceptibility, x(/i) at a renor- 
malization scale (/x = IGeV^) in the following manner: 

cl){u) = x{^)u{l-u). (39) 
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The remaining functions ip^'^\u), ip^^\u), A{u), and B(m) are also defined as |l2l. Il3|: 



15 



ifj'' \u) = -20m(1 - u){2u - 1) + ^{u^ - 3u}^)u{l - u){2u - 1)(7(2m - 1)^ - 3), 



16 



19 



tlj'-'iu) = (1 - {2u - ir){5{2u - ir - 1)^(1 + ^u}; 



16^-^' 



A(u) = 40m(1 - u){3k - A;+ + 1) + 8(^2^ - 8^2) 

X [m(1 - m)(2 + 13m(1 - u)) + 2m3(10 - 15m + 16m=^) \nu 



+2(1 - uf{10 - 15(1 -u) + 6(1 - M^)) ln(l - u)], 



B(u) 



40 



da{A - a){l + 3A;+)[-- + -{2a - 1] 



(40) 



where k, ^2, ^2" /s^ are constants (see [12|, [ISJ). Putting the above equations all 
together after performing the integrals over d'^x and d'^k, the coefficients of the corresponding 
structures ie^yape°'q^ and [q^e^ — e^q^] are obtained as follows: 

nnon—pertf \ Qs f 

If {P,q) = — ^ / 



(i'u|mj?(ss) A(m) 
mc{ss)B{u) - 2(5m^ - p'^){mc{ss)4>{u) - f^^i)^^\u)) | 



-r-rnora— pert 



(iu<^ A(M)m^(ss) + 2(-5m^ +p^)(ss) 0(m) 



2(m2 — p2)3 



u 



(41) 



Now, To find the QCD sum rules for the form factors, we match the coefficient of the 
selected structures from both phenomenological and QCD sides and perform the Borel trans- 
formation with respect to the momentum of Dsi meson, (p2 Ml). To further suppress the 
contribution of the higher states and continuum, we also perform the continuum subtraction 
and use the quark-hadron duality assumption. As a result, we get 



jDsi J(ms- 



{nis+mc) 



, S-Q2 + ^^1,: 



,2(c+d+e+f) 



(42) 



where Sq is the continuum threshold; and V and A correspond to 1 and 2 in the right hand 
side, respectively. To obtain the expression of the above sum rules in Borel scheme, we 
perform the Borel transformation using the following standard rule: 



B 



(p2 - sy 



r(n)(M2)'^-i' 



(43) 
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IV. QCD SUM RULES FOR THE FORM FACTORS RESPONSIBLE FOR THE 
ELECTROMAGNETIC PENGUIN MODE 



At quark level, the FCNC based EP transition of Be — > Dgi'f proceeds via b — > 57, whose 
effective Hamiltonian is written as 



Heff 



Gpe 



The amplitude of this mode is obtained from 



1 + 75 , 1-75 
nib — 7: 1" TT^s 



hence to proceed further, we need to calculate the following matrix elements: 



(L',l|s(7^,(l±75)0|5c), 



(44) 



(45) 



(46) 



which can be parametrized in terms of two gauge invariant form factors Ti(g^) and T2{q^) 
in the case of real photon, i.e. 



{Dsi{p,s^''^''^)\sa,,q''jrMBcm = ^WAe^""'^^ VQ'Ti(O), 



(47) 



These two from factors are not independent of each other and using the relation afj,,^^^ — 
— ^^nuapo-"'^ we read Ti(0) = ^T2{0). Therefor, we need to calculated just one of them and 
in the following, we calculate the 22(0). The corresponding correlation function is chosen as 

^;..(p^g') =^'y" j ciWye-(«--^'^)(0|r{c(t/)7.75s(y)s(0)(7^„g"6(0)6(x)75c(a;)}|0), 

(48) 

where, 675c and c'jujss are the interpolating currents of the initial and final mesonic states, 
respectively; and safj^aq'^b is the transition current. Prom the general philosophy of the QCD 
sum rules, we calculate this correlation function again in two different languages, namely 
hadronic and quark-gluon. Por the hadronic or phenomenological part, we get 



{m%^ - (52) {m%^^ - p2) (rrib + rric) 



{mB^-mjjJgi,^T2{0) 



m 



)p^p.T2{0) + {p + Q)^ 



p.q 
—Pv - qy 



m 



Ds 



r2(o)U...., (49) 
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where ... denotes contribution of the higher states and continuum which will be suppressed 
then applying Borel transformation as well as continuum subtraction. To calculate the 72(0), 
we choose the structure g^j,^. In deriving the above equation, we have used the following 
definition of the decay constant of the Be meson: 

{Bc\Hc\0) = ^T^^x^- (^0) 
{mb + rric) 

In QCD side, the correlation function is written in terms of the selected structure as 

^^. = 5;..^(p^g2), (51) 

where 

Here the perturbative part is related to the spectral density, p^^^* by the double dispersion 
integral, 

Q^) = j I d^'dtj£^^0^^ + subtraction terms, (53) 

and for the non-perturbative contributions, we will calculate the two-gluon condensate dia- 
grams. 

Now, we focus our attention to calculate the spectral density. Using Cutkosky method 
to calculate the spectral density, we get: 

pP^'^(s', t) = 2A^c|/o A[{mb - mc){mc + m^) + - A'[(m6 - ■mc){mc + m^) + s'] 

+2mc[{mc + ms)s' + m^t — rrict] — mc{mb + ms)u + 2^41 (—2s' + m)|, (54) 

where 



A'(a, 6, c) = + 6^ + - 2ab - 2ac - 26c, 

A, = , ~/l ,J AH + A'h' - AA'u + mli-As't + u% 

{—As't + w^Y 

A = s' + ml — ml, 
A' ^t + ml-ml, 
u^t + s' -q'^. (55) 
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FIG. 3. Feynman diagrams for gluon condensates corrections. 

The following inequalities, that comes from the (5-integration of the Cutkosky method, re- 
strict the boundaries of the s'- and i-integrals: 

ml <t < to, 

t-J^<s'<s',. (56) 

There are several sources for non-perturbative contributions, such as quark, quark-gluon, 
and gluon condensates. But the quark condensates (s-quark) and quark-gluon condensates 
give zero contributions after applying the double Borel transformation with respect to 
{Q"^ — > Mf ) and {p^ — > M|). Therefore, the remaining source of non-perturbative contri- 
butions would be the gluon condensates [see Fig (3)]. The calculation of such contribution 
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is lengthy but standard. For the nonperturbative part in Borel scheme, we get 



^non-pert ^ ^2^2 /^q2\(J^^^ (57) 



TC 



where Cq^ is the Wilson coefficient of the gluon condensates that is defined as 

Cq2 = Cq2 + Cq2 + Cq2 + Cq2 + Cq2 + C*^2- (^S) 

The explicit expression of the is given in the Appendix. 

From the same procedure as presented in the previous section, we get the following 
expression for the form factor under consideration: 

/ I \ ml /M? m?n l^j 



7^2(0) 



(2n) 




(59) 



where Sq and to are the continuum thresholds in the initial and final channels. 

V. NUMERICAL ANALYSIS 

This section is devoted to the numerical analysis of the form factors as well as estimating 
the branching ratio in WA and EP channels as well as the total branching fraction. For this 
aim, we use the quark and mesons' masses as rric = (1.275±0.015) GeV , rris ^ 142 MeV jl^ . 
mb = (4.7±0.1) GeV 15|, m^.i = (2459.6±0.6) MeV, niB, = (6.277±0.006) GeV Q. For 



the va^ 
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ues of the decay constants, we use fn , = (225 ±25) MeV and /^^ = (350 ±25) MeV 
19|. The values of the condensates are [jj: {ip^p\f,=iGev) = -(240 ± 10Me\/)^ {ss) = 



(0.8 ± 0.2) (^V^), ml = (0.8 ± 0.2) GeV^ and {^G^) = 0.012. The parameters entered the 
photon DAs are also taken as x = 3.15 ± O.SGeV-^, k = 0.2, fc+ = 0, Ci = 0.4 Ci_= 0, 
C2 = 0.3, (2^ = 0, = -(4±2) X IQ-^GeV^, = -2.1±1.0 and u}^ = 3.8±1.8 QQEl- 
The remaining parameters are chosen as \Vcs\ = 0.957±0.017, \Vcb\ = 0.0416±0.0006, \Vtb\ = 
0-77l|J:24, \Vts\ = (40.6±2.7) x 10-^ Gr{fx = m,) = -0.0068-0.02i, and r^, = 0.52 x lO'^^^. 

The sum rules for the form factors contain also continuum thresholds and Borel mass 
parameters as auxiliary objects. We should find their working region such that the physical 
observables be practically independent of them. The continuum thresholds are not totally 
arbitrary but they are correlated with the energy of the first excited state in the initial and 
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final mesonic channels. Our numerical results show that, the results depend weakly on the 
thresholds in the intervals, sq = to = {6 — 8)GeV'^ and Sq = (45 — 50)GeV'^. The working 
regions for the Borel parameters are obtained demanding that not only the contribution 
of the higher states and continuum are effectively suppressed but also the contribution of 
the higher order operators and higher twist DAs are small, i.e., series of the sum rules are 
convergent. These conditions lead to the intervals 6 GeV^ < < 12 GeV^, 10 GeV^ < 
Ml < 30 GeV^ and 5 GeV^ < M| < 12 GeV^ for the Borel mass parameters. 

Now, we proceed to find the fit functions of the form factors using the aforesaid working 
regions for the auxiliary parameters as well as other input parameters. Here, we would like 
to mention that for the decay rates we need the values of the form factors Fy^"^^ and F^"^^ 
at = m\^^ for Flf"^ and F^"^ at = rnj^^^ and for T2 at = 0. However we determine 
their fit functions in general then give the their values at fixed points. The fit functions for 
the form factors Fy^"^^ and F^^"^^ are given as 

/(g2) = J(^) (60) 

"si "si 

where /(O), a and h are the fit parameters whose values are given as: 



form factors 


/(O) 


a 


h 


4^-)(Q2) 
Ff-)(Q2) 


0.134 

-1.355 


0.145 
3.217 


-0.008 
-0.205 



The values of these form factors at = m^^ are obtained as: 



m%J = 0.083 



ADsi) 



m 



Be 



-0.102 



Also the fit functions for the form factors F^y are as follows j4| 

Fv{0) 



V 



1 — /m\ 

Fa{0) 
1 — p'^/ml 



(61) 



(62) 
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where the fit parameters are given by 



Fy(0) = 0.44GeV 


ml = 43.1GeV^ 


F^(0) = 0.21GeV 


ml = 48.0GeV2 



The values of the form factors -Fy^"* obtained at = "^d^i given as: 



V 



m 



0.512GeV 



Fj^^) (p2 ^ ^2^^^ ) ^ 0.240GeV (63) 

For the form factor induced by the electromagnetic penguin at = 0, we obtain 

7^2(0) = -0.298 (64) 

At the end of this section, we would like to calculate the total decay width and branching 
ratios. Using the amplitudes of each decay mode, we find the following expressions for the 
total decay rates at fixed points in different channels: 
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Be 



ml 



m. 



ml 



(65) 



m. 



10247r4 

X (^17m^ + SOmbm, + 17m^) [T2(0)]^ 



Gla /ml 



mf) \3 



647r^|K,K*/ 



10247r4 



(66) 



+ 



Be 



Be 



y 

m'^ 



Bo 



Mc^Y\Vtjyi\\nmi - mm^ms + nmi)[T2 



+i67r^r2(o)|K6v;:ilW,| 



/D,ir«i:,,iX{4 — —{mb - m^) 



mt 



(Be) 
V 



m 



-(mb + ms)} + /B,{Fr°'^(m6 + m,)X 



B. 

(Dsi) 



4Fp'\mb - ms)Y} 



(67) 
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where X and y, are the real and imaginary parts of the Wilson coefficient C7, respectively. 
In these formulas, the fixed point values of the form factors discussed in the above are used. 

Finally the numerical values of the corresponding branching ratios for the radiative decay 
under consideration are obtained as follows: 

B^^^\B^ D,a) = 1-765 X 10-^ 
B^^^\B^ Dsi-i) = 3.034 X 10-^ 

g(totaO(5^ ^ D^^^) = 3.136 X 10"^ (68) 

The order of the total branching ratio indicates that this decay channel can be detected at 
LHC. Any measurement in this respect and the comparison of the obtained data with our 
predictions in the present work can give valuable information about the nature and internal 
structure of the participating particles especially the Dgi meson. 
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VI. APPENDIX 

The explicit expressions for C^2 are given as follows: 

eg, = m,{2ml{Io[l, 3, 1] - 3mllo[l, 4, 1] + 2m^/o[l, 5, 1]) 

-2mbml[Io[l, 3, 1] + {-3ml + ml)Io[l, 4, 1] + 2m^(mb - m,)(m6 + m,)/o[l, 5, 1]] 
+2me(/o[l, 2, 1] - 4m^7o[l, 3, 1] - m^m./ofl, 3, 1] + m^/ofl, 3, 1] + 5m^7o[l, 4, 1] 

+3mlmJo[l,A,l]-3mlm%[l,A,l]-mbm%[l,A,l] 

+2ml{-rrn, + m,)(mfe + m,)2/o[l, 5, 1]) + 2m,[/o[l, 2, 1] - 2m^/o[l, 5, 1] + ij^'^^l, 3, 1] 
+mt{5Io[l, 4, 1] + 2inh 5, 1]) + m^(-47o[l, 3, 1] - 34'^'^1[1, 4, 1] + 3t'^[l, 3, 1])] 
+m,m^[/f ^1[1,4, 1] + 3/f '°][1,4, 1] - 2m^(/;^'^] [1, 5, 1] + 3/f '°1[1,5, 1])] 
+m,(3/^^5 [1, 3, 1] + J^'l [1, 4, 1] + 4^'°] [1, 3, 1] + 2m^(3/^^l [1, 5, 1] + 4^'°^ [1, 5, 1]) 
-mUgin^, 4, 1] + 24°'^1[1, 5, 1] + 34^'°] [1, 4, 1] - 24^'°1[1, 4, 1]) 
-4/f^l[l,5,l] + 44^'°][l,3,l])}, (69) 
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-6mbm1lo[l, 1, 1] + 6mcm^/o[l, 1, 1] - 5m^m^/o[l, 1, 2] - 8mbmcm1lo[l, 1, 2] 

+5mlmcm%[l, 1, 2] + m>2jp[i^ 2] + 2mf,m^/o[l, 1, 2] - 2mem;Vo[l, 1, 2] 
-12mbmcm%[l, 1, 2] + 12m>^/g[i^ 1^ 2] + llm6mj/o[l, 1, 2] - llmem^/o[l, 1, 2] 
-7mlmcmllo[l, 1, 3] + lmlmlmlh\\, 1, 3] + m^m^7o[l, 1, 3] - mlmjn%{l, 1, 3] 
-10m^mcm^/o[l, 1, 3] + 10mgm^m^/o[l, 1, 3] + 6m^m^Jo[l, 1, 3] 
+1 7mbmcm^ /o[l, 1, 3] — 6m^mcm^/o[l, 1, 3] — 17m^mj/o[l, 1, 3] + mbm^/o[l, 1, 3] 

mlhW 1, 3] - 22m>^7o[l, 1, 3] + 6mem%[l, 1, 3] 
-2mbmJ/o[l,l,3] +m3jo[l,l,4] +me/o[l,l,4] +mfeme/o[l,l,4] -mgm,/o[l,l,4] 
-m^ms/o[l, 1, 4] - mjnshW, 1, 4] + 2mbmcmsIo[^, 1, 4] + m^mcms/o[l, 1, 4] 
+10m^mem^7o[l, 1, 4] - 10mlmlm%[l, 1, 4] + 2mlm%[l, 1, 4] - 2mlmcm%[l, 1, 4] 
+12mfmcm^/o[l, 1,4] - 12mgm^mf/o[l, 1,4] - 10m6mcm^/o[l, 1,4] - Unibm ^ml I o[l, 1,4] 
+12m^mJ/o[l, 1, 4] - mjoll, 1, 5] - mbmJo[l, 1, 5] + mlmJo[l, 1, 5] - mlm%[l, 1, 5] 
+mcms/o[l, 1, 5] - 2mbmcmslQ[l, 1, 5] + 2mlmcmslQ[l, 1, 5] - 2mlmlmsh[l, 1, 5] 
-6ml/^^] [1, 1,2]- m,J^^] [1,1,4]+ m jf^' [1,1,4] + 21/2m^/f ^1 [1,1,4] 
+bmbmllf^\l, 1,4] - m,m^/^^l[l, 1,4] - 13/2m^/^^l[l, 1,4] - 3/2/J''^][l, 1,5] 
+3/2m24°'^l[l, 1, 5] + m,m4°'^l[l, 1, 5] - 21/2m^m2/f '^^[1, 1, 5] + ISm^^'^lfl, 1, 5] 
-15m^m^^'^I[l, 1, 5] - 6m,m^/^^l[l, 1, 5] + 15m^m^/^^l[l, 1, 5] - 7m^^'^l[l, 1, 5] 
-2m,m^/£''^] [1, 1, 5] + ISmgm^/^^l [1, 1, 5] - nis^'^ [1,1,2]- 6m^/^'l [1, 1, 2] 
+7/2m24°'2l[l, 1, 3] + 5m34°'2l[l, 1, 3] + 7/2m^4^'°l[l, 1, 2] + 3/24^'°][l, 1, 3] 
+m67f 1, 3] + 3mjf '"Jfl, 1, 3] - Smfcm^/f 1, 3] - 7/2m^m^7f 1, 3] 
+m,m^4^'°l [1, 1, 3] - 5mt4^'°l [1, 1, 3] - Qml^'^ [1,1,3] + QmlmlC^ [1, 1, 3] 
-18ml4^'°][l, 1, 3] + l/2m^4^'°][l, 1, 4] - m,m,4^'°l[l, 1, 4] + 5m^4^'°l[l, 1, 4] 
+51/2mt7f 1, 4] + 6m,mt7f 1, 4] + 33m^4^'°][l, 1, 4] - 5m^7f '"^[l, 1, 3] 
+l/27r^[l, 1,4] + m,7r^[l, 1,4] + 6m,^7r^[l, 1,4], (70) 
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l/6{2m^m27o[3, 1, 2] + 2m^m,(7o[3, 1, 1] + m^7o[3, 1, 2] - 7o[3, 2, 1]) 

+2ml[-mcm%[3, 1, 2] + m^/o[3, 1, 2] + m,(-/o[3, 1, 1] + m^/o[3, 1, 2] + ij^'^^S, 1, 1]) 
+m,K/o[3, 2, 2] + /^''^][3, 1, 2])] - 3/^^][3, 2, 1] + 3m^/^'^][3, 2, 1] + /^^'^'[S, 2, 2] 
-2m^4°'^l[3, 2, 2] + m^7o * [0, 1][3, 2, 2] + 3/^^(3, 1, 2] - 2/^^(3, 2, 1] + 2m^/r^[3, 2, 2] 

[3, 2, 2] + 3jf [3, 1,1]- 6m^/f [3, 1, 1] - ij^''^ [3, 1, 2] 
+2m3m,{/o[3, 2, 1] + m%[3, 2, 1] - /o[3, 2, 2] - 4/^^'^l[3, 1, 2] + 3/J'''][3, 2, 1] + 3/f '°][3, 1, 1] 
-21^3, 2, 1] + m^[-27o[3, 1, 2] + 7o[3, 2, 2] - 3(7^ ■^1[3, 2, 2] + 7^'°1[3, 2, 1])]} 
+3m^7f '°1[3, 2, 2] - 7^ '^'[3, 1,1] + 7^ '^'[3, 1, 2] - m^7f '^^[3, 1, 2] 
-2m3[m^7o[3, 2, 1] - m,mtlo[3, 2, 2] - 2m3(7o[3, 1, 1] - 7^ '^1[3, 2, 1] + 7f '"l [3, 2, 1]) 
+mem^(2m^7o[3, 2, 2] - 1^^, 2, 2] + 37^ '°][3, 2, 2]) - me(m^7o[3, 1, 2] - 27^^'^J[3, 2, 2] 
+7f '°1[3, 1, 2] + m^(7o[3, 1, 2] - 27o[3, 2, 2] - 3(7^''^1[3, 2, 2] + 7^ '°][3, 2, 2])) + 7^ '^1[3, 2, 1]) 
+m,((l + 2m^)7o[3, 2, 2] + 27^^! [3, 1, 1] - 27^''^] [3, 2, 1] + Amlli^^'\3, 2, 2] + li^^'\3, 1, 2] 
+2mlin3, 1, 1] + 27i^'°l[3, 2, 2] - 27^ '^][3, 2, 2])] + 2m^7f ■^1[3, 2, 2] + m^(2m^7o[3, 1, 2] 
+2mX^[3, 2, 2] + 2m,m,(m^7o[3, 1, 2] - 7o[3, 2, 2] + 1^^, 1, 1]) - 37^^''][3, 2, 2] 
+37^''^] [3, 2, 2] - 7f '°][3, 1, 2] + m^(-2m^7o[3, 1, 1] - 37^'^! [3, 2, 2] + 7f '°][3, 1, 2]) 
+m^(37^^'^][3, 1, 2] + 1^3, 2, 2]) + 7f '^^[3, 2, 2]) + 2mem,[-7^''^][3, 2, 1] 
+7f '^][3, 2, 2] + 7f ■°1[3, 2, 2] - 7f ■^1[3, 1, 2] + m^(7f '^][3, 2, 2] - 27? ■°1[3, 1, 1] + 7? ■^1[3, 2, 2])] 
-7? '^'[3, 2, 2] - 7f '^Ip, 1, 2] + 27f °][3, 1, 2] + 6m^7f °][3, 1, 2] - m^(2mt7o[3, 1, 1] 
+24°'^1[3, 2, 2] + 4°''1[3, 2, 2] - 91^^, 1, 2] + 77? '°][3, 2, 2] + m2(27o[3, 1, 2] 
-27o[3, 2, 1] + 77? '^][3, 2, 2] + 97?'°1[3, 2, 2]) - 7?'°1[3, 2, 2]) - 37?'°][3, 2, 2] - 2m^4''°l[3, 2, 2] 
-3m^7?'°'[3, 2, 2] - 2mb{mlmllo[3, 2, 2] + m>,(m^7o[3, 2, 2] + 7?''' [3, 1, 2]) 
-ml[{l + 2ml)Io[3, 1, 2] - 7o[3, 2, 2] + 37?'^][3, 1, 2] - 27?'^][3, 2, 2] - 47?'°1[3, 1, 2] 
+m^(-7o[3, 2, 1] + 7?'^1[3, 2, 1] + 57?'°][3, 2, 1]) + 37?'°1[3, 2, 2]] + m>,((l + m^)7o[3, 1, 2] 
-(1 + 2m^)7o [3, 2, 1] - 47?'^^ [3, 1,2] + 27?'^1 [3, 2, 1] - 2m^7?'^l [3, 2, 1] - 37?'^^ [3, 2, 2] 
-4m^7?'°l[3, 1, 2] + 3^'°I[3, 2, 1]) + m,(-7?'^][3, 2, 2] + m^(7o[3, 2, 1] - 7?''^][3, 2, 1]) 
+7?'°][3, 2, 2] + 7?'^1[3, 1, 2] - 7?'°1[3, 2, 2] + m^(7o[3, 1, 1] - 7o[3, 2, 1] + 7?'^][3, 2, 1] 
-37?'°1 [3, 1, 2] + 7?'°] [3, 2, 1] - 7?'^] [3, 1, 2] + 27?'°] [3, 2, 2])) - m,(7?''l [3, 2, 2] 
+27?'°] [3, 1, 2] + m^(7o[3, 1, 2] - 7?'°][3, 2, 2]) - 7?'°][3, 2, 2] - 7?'°][3, 1, 2] (71) 
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+m^(7o[3, 2, 1] - 27o[3, 2, 2] + 27^ '^1[3, 1, 2] + 27? '°1[3, 1, 1] - 27? '^][3, 2, 2] + 7? '°1[3, 2, 1]) 

+7?'^] [3, 1, 2])) - 3m27f ^][3, 2, 1] - m^(2m^7o[3, 2, 2] + 7?'^1[3, 1, 2] - 67?'^' [3, 2, 1] 
-27?'^][3, 1, 2] + 64°'^1[3, 2, 2] + 7?'^1[3, 2, 2] + 3m^7f '"llS, 1, 2] + 2m>,(7o[3, 1, 2] 
+2(-l + m^)7o[3, 2, 1] - 57^^'^1[3, 2, 2] - 7? '°1[3, 2, 1]) + 37? '°1[3, 2, 1] - 27? '°1[3, 2, 2] 
-7m^7?'°l[3, 2, 2] + 2m,m,(7o[3, 1, 2] + (-1 + m^)7o[3, 2, 1] - 2m27o[3, 2, 2] 
+m^7o[3, 2, 2] + 3m^7?'^][3, 1, 2] + 27?'^][3, 1, 2] - m^7?'°][3, 2, 1] - 7?'^^ [3, 2, 1]) + 27?'^1 [3, 2, 1] 
-37?'^1[3, 2, 2] + 27?'°1[3, 2, 1] - 2m^(7?'^l[3, 1, 1] - 7?'^1[3, 2, 2] - 7?'°1[3, 1, 2] + 37?'°][3, 2, 1] 
+7?'°1[3, 2, 1]) - 2m^(m27o[3, 2, 2] + m^7o[3, 2, 2] + 27^''^1[3, 2, 2] - 57ri[3, 2, 2] - 27?'°][3, 2, 2] 
-77?'^] [3, 2, 1] - 7?'°] [3, 2, 2]) - 27?'^1[3, 2, 2]) + 7f^l[3, 2, 2] + m^(-27;^'^^[3, 2, 2] 
+107?'^][3, 1, 2] - 27?'^][3, 2, 1] + 7?'°1[3, 2, 1] - m^(7?'^l[3, 1, 1] - 37?'°][3, 2, 2]) + 7?'°][3, 2, 2] 
-147?''1[3, 1, 2] + 147?''' [3, 2, 2] + 37?''^ [3, 2, 2] + 27?'°1[3, 1, 2] - 107?'°' [3, 2, 1] 
+2m^(7?'^l[3, 1, 1] + 7?'^1[3, 1, 2] + 7?'^][3, 1, 2] - 27?'°1[3, 1, 2] + 77?'^][3, 2, 2] + 57?'°1[3, 2, 2]) 
-37?'^1[3, 2, 2] - 37ri[3, 2, 2]) - 7?''^1[3, 2, 2] + m^7?'°][3, 2, 2] + 7?'^1[3, 2, 2]}, 

= l/12{2m^((-l + m2)7o[3, 1, 2] + 7o[3, 2, 1]) + 2m^m,(7o[3, 1, 1] + m^7o[3, 2, 1]) 
+24m^K + ^s)/o[3, 2, 2] + 7?''' [3, 2, 2] + 2m^(8mem,7o[3, 1, 1] + 8m27o[3, 2, 2] 
+37?''1[3, 1, 2] + 7?'°1[3, 1, 1]) + 6m^(2m^7o[3, 1, 1] + 2m>Jo[3, 2, 2] - 2m,(-27o[3, 1, 2] 
+(5 + 2m^)7o[3, 2, 1] + 27?''^ [3, 2, 2]) + me(-2(3 + 2m^)7o[3, 2, 2] - 97?''^ [3, 2, 1] 
+7?'°][3, 1, 2])) + 37?'°1[3, 2, 1] - 3m27?'°l[3, 2, 1] - 7?'°1[3, 2, 2] - 2m,m,(27?'^l[3, 1, 2] 
-7?'°1 [3, 1, 1] - 7?'°][3, 1, 2] + m^(7o[3, 2, 2] + 7?'°][3, 2, 2])) + m^7?'^l[3, 1, 2] 
+27?'^1[3, 2, 2] + 7?''1[3, 1, 2] + 3m^7?'°l[3, 1, 2] - m^(2(-l + m^)7o[3, 1, 2] + 27o[3, 2, 2] 
+m^7i"''] [3, 1,1] + 27i^''l [3, 2,2] + 1^'''^ [3, 2,2]- 57?'°1 [3, 2, 1] + 37?'°' [3, 2, 2] 
+5m27?'°l[3, 2, 2] - 7?'°1[3, 2, 2]) - 7ri[3, 2, 2] + m^(4m^7o[3, 2, 1] + 4m>,7o[3, 2, 2] 
+24°'^1[3, 2, 1] - 157?''1[3, 2, 2] - 27?''][3, 1, 1] - 4mem,(4(-l + m2)7o[3, 2, 1] 
+57^^'^] [3, 2, 2] - 47?'°1 [3,1,2]) + 67?'°^ [3, 2, 2] - 6m^7?'°l [3, 2, 2] 

-2m^((6 + 8m^)7o[3, 2, 2] + 3(7j''^][3, 2, 2] + 7?'°][3, 2, 2])) + 67?'^1[3, 1, 2] + 47?'°1[3, 2, 2]) 

+3m6(2m>,(-7o[3, 2, 1] + m27o[3, 2, 2] - 27?'^][3, 2, 2] + 27?'°1[3, 2, 2]) 

+me(-4(-l + m^)7o[3, 1, 2] - 47o[3, 2, 2] - 3m^7?''][3, 1, 1] + 37?''1[3, 2, 1] (72) 
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-9irV, 2, 2] - 37r^[3, 1, 2] + 57i''°J[3, 1, 2] + 1^^, 2, 1] - Bmpn^, 2, 2] 

+3/f °1[3, 1,2])- 4m,(/o[3, 1, 2] + (-1 + m2)/o[3, 2, 2] + /^'''1[3, 2, 2] 
+/f '°][3, 1, 2] - 24^'°][3, 2, 1] - m^/f '°][3, 2, 2] - /f '^^[3, 2, 1] + I^^''\3, 2, 2])) 
-3mf (4m^(m^/o[3, 1, 1] - /o[3, 2, 2]) + 2m>,((3 + 2m^)7o[3, 1, 2] - il^^, 2, 2] 
+4/f [3, 1, 2]) - 4m,(-3/o[3, 2, 1] + 4Jo[3, 2, 2] + 3jJ'''l [3, 2, 2] + m2(3/o[3, 1, 2] 
-2/f '°][3, 2, 2]) - 44^'°] [3, 2, 2] - 2/f '^'[3, 2, 2] + 2/^^'°][3, 2, 2]) + m,(-16/o[3, 1, 2] 
+127o[3, 2, 2] - 277^''] [3, 1, 2] + 67i°''l[3, 2, 1] - 64°'"" [3, 2, 1] + 31^^, 1> 1] 
-2m^(6/o[3,2,l] + 3/^^l[3,2,l] + 54^'°1[3,1,2]) + lo4^'°][3, 2, 2] + 64^'°1[3, 2, 2])) 
+ml{mt{-6Io[S, 1, 2] - 4m2/o[3, 2, 1] + 4Jo[3, 2, 2]) + m3(4m,/o[3, 1, 2] - 6m,/o[3, 2, 1] 
-4mf7o[3, 2, 2]) - 37^ '^][3, 1, 1] + 12/f '^^[3, 1, 2] - 91^^, 1, 2] + 47? '°][3, 2, 1] 
+2m,m,(77o[3, 1, 1] + 97^'''] [3, 2, 2] + 2m2(37o[3, 1, 2] + 7? •°1[3, 1, 2]) - 27i''°][3, 2, 1] 
-67f '°][3, 2, 2]) + 9m,27f '°][3, 2, 2] - 97?'^^ [3, 1, 1] - 2771',!^^ 2, 2] - 27?'^^ [3, 2, 1] 
+m^(147o[3, 1, 2] - 127o[3, 2, 2] + 97^^''1[3, 1, 2] - 27? ''][3, 2, 2] + 27?''1[3, 1, 2] 
-107?'°1[3, 1, 2] + 2m^(67o[3, 2, 1] + 7?' ^'[3, 2, 1] + 57? '°1[3, 1, 2]) + 97?'°^ [3, 2, 1] 
-27?'°1 [3,1,2])- 67?'°1 [3, 1,2] + 67?'°^ [3, 2, 2] - 6m^7?'°] [3, 2,2] + 24^'°1 [3, 1,0]) 
-4''°'[3,2,2]}, 



= l/6{-2mcm,(m^(27o[l, 2, 2] - 37o[l, 3, 1]) - 7o[l, 3, 2] + 2m^7o[l, 3, 2] + 7o[l, 3, 3]) 
+4m^(2m,m27o[l, 2, 1] + 4m37o[l, 2, 1] - mJo[l, 2, 3] - 2m,7o[l, 3, 3]) 
-2ml[{l + 2m2)7o[l, 2, 2] - 3m,27o[l, 3, 1] - 7o[l, 3, 3] + 2m^7o[l, 3, 3]] - 37?''! [1, 1, 3] 
+7?'^][1, 2, 2] + 6m^7?'^][l, 2, 3] + 2m^7?'^][l, 2, 3] - 3m^7?'^][l, 3, 2] + 2m^7r][l, 1, 2] 
-7?'^1 [1, 2, 3] - 9m^7?'°] [1, 1, 3] - 7?'°] [1, 2, 2] + 6m^7?'°] [1, 2, 2] + 37?'°] [1, 3, 1] 
-2m^(me(-2m,7o[l, 2, 1] + 4m^7o[l, 2, 2]) + ml{4:m%[l, 2, 1] - 27o[l, 3, 3]) 
+(3 - 6m^)7?'^][l, 3, 3] + 7?'°1[1, 2, 2] - 2m^7?'°][l, 3, 2]) - 2m^K(-37o[l, 2, 3] 
-6m2(7o[l, 1, 2] - 7o[l, 3, 2]) + 4m,i7o[l, 3, 2] + 27o[l, 3, 3]) + 2m,(4rr47o[l, 2, 3] 
-7o[l, 3, 3] - 27?''1[1, 2, 2] + m2(-67o[l, 1, 2] + 67o[l, 3, 2] + 47?''][1, 3, 3] - 47?'°1[1, 1, 3]) 
+27?'°1[1, 3, 2])] + 2m,{m,[{-l + 2m2)7o[l, 2, 2] + 2m^7o[l, 2, 3] - 3m^7o[l, 3, 1] 
+7o[l, 3, 2]] + 2m,[2m^7o[l, 2, 2] - 7o[l, 2, 3] + 7o[l, 3, 2] - 7?''1[1, 1, 3] + 7?'°1[1, 2, 3] (73) 

23 



+m^(27o[l, 2, 2] - 37o[l, 3, 1] + 21^^, 2, 3] - 21^^, 3, 3])]} + 2mpnh 3, 3] 

-2m2/f 1,3] + 3, 2] + ml{2ml{4m%[l, 1, 3] + 2/o[l, 2, 3] - 6m2(/o[l, 1, 2] 

-/o[l, 3, 2]) - 3/o[l, 3, 2]) + 2m,ms{2Io[l, 1, 2] - 6ml{Io[l, 1, 3] - /o[l, 3, 2]) - 3/o[l, 3, 3] 
+Am%[l, 3, 3]) + (9 + 6ml)I^^''\l, 2, 3] - 27^^'^1[1, 3, 2] - 18m'/^'% 3, 3] - 4m^^'^l[l, 3, 3] 
-4m^/[°'2l[l, 2, 1] + 2^'^][1, 3, 3] - 61^1, 1,3]- 6m^4^'°][l, 2, 2] + 18m^/f '"'[l, 2, 3] 
-Umtinh 2, 3] + 3/f 3, 2] + 4m^/f '"'[I, 2, 2] - 2/f '°1[1, 3, 3])}, 



C^, = 2/3m,{m3[2m^/o[2, 1, 3] - (-1 + ml)Io[2, 1, 4] + m2{-3/o[2, 1, 2] + 2(m2/o[2, 1, 3] 
+^ '^^[2, 1, 4])}] + mlm,[-2m%[2, 1, 3] + (-1 + m^)/o[2, 1, 4] - 21^^, 1, 4] 
+m^(37o[2, 1, 2] - 2m^7o[2, 1, 3] + 47^^''1[2, 1, 4] + 2/f '°][2, 1, 4])] + mf,[-2mpo[2, 1, 4] 
+1^2, 1,3] - 24°'^1[2, 1,4] + in^, 1,2] + /r'[2, 1, 3] + m^/f '°1[2, 1,3] 
+m^(5/o[2, 1, 2] + 2/f [2, 1, 3]) - /f '^1 [2, 1, 4] + m^(-4/o[2, 1, 1] - 2/f^l [2, 1, 3] 
+64°'^1[2, 1, 4] + m^(37o[2, 1, 2] + 27^^'^^ [2, 1, 4]) - 34^'°1[2, 1, 4] + 27^'^] [2, 1, 4])] 
+me[27f ^1 [2, 1, 3] + 7^'''] [2, 1, 3] + m^(-57o[2, 1, 2] + 2m27o[2, 1, 3] - 67^'°^ [2, 1, 4]) 
-37f [2, 1, 4] - 7f °] [2, 1, 4] + m^(47o[2, 1, 1] - 3m^7o[2, 1, 2] - 47f''^l [2, 1, 4] 
-27r^[2, 1, 2] + 91^2, 1, 3] + 2^ ■°1[2, 1, 4])]}, (74) 

where 7„[a, b, c] and 7n'''^[a, b, c] are defined as: 

(_-]\a+b+c 

''^t"' ^' ^] = 167r^r(a)r(b)r(c) ^^'^"""'^^^'^'"""'^ t/o(a + 6 + c - 4, 1 - c - 6), 

|'_-|\a+b+c+l 

^■I"' '1 = i6.T(;)r(6)r(c) '^-)"°"'^-'°'°" t/„(a + . + c - 5, 1 - c - 6), 

'1 = 16.^"r(a)r(6)r(c) '^^?>"°"<^'^^^'"°" C/„(a + + c - 5, l - c - 6), 
7M[a, 6, c] = ^ 4- (75) 

The function f/o(fl, b) in the above formulas is given by 

Uo{a,b) = [\y{y + M^ + M^ry'eM-^-Bo-B,y], (76) 

Jo y 
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where 



Bo 



m 



[Mlml + Ml{ml + ml)i 



(77) 
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